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I. Introduction 

The asymmetric simple exclusion process (ASEP) is a special case of processes 
introduced in 1970 by F. Spitzer jS]. In ASEP, particles are at integer sites on the 
line. Each particle waits exponential time, and then 



(1) with probability p it moves one step to the right if the site is unoccupied, 

CN ■ otherwise it does not move; 
> 

CN ! (2) with probability q — 1 — p it moves one step to the left if the site is unoccupied, 

otherwise it does not move. 

In the totally asymmetric simple exclusion process (TASEP) particles can only 
move in one direction, so either p = or q = 0. In a major breakthrough, K. Johans- 
son [2] related a probability in TASEP to a probability in random matrix theory. If 
q = and initially particles were at the negative integers, then the probability that 
at time t the particle initially at — m has moved at least n (> m) times equals the 
probability distribution for the largest eigenvalue in the Laguerre ensemble of m x m 
matrices with weight function x n ~ m e~ x . Thus, it is given by a constant depending 
on m and n times the determinant 

det ( I x n - m +*+i e ~ x dx 



i,j=0,...,m-l 

This connection led to considerable progress in understanding TASEP and the deriva- 
tion of asymptotic results. For ASEP there is no longer a determinantal structure 
and a different approach was required. 



1 This is an expanded version of a series of lectures delivered by the second author at Universite 
de Paris in June, 2009, describing the results in the articles [9l HOI HTJ . Although complete proofs 
will in general not be presented here, at least the main elements of them will be. 
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Here are the main results of [91 HOj [TT] . First we consider ASEP with finitely many 
particles. For iV-particle ASEP a possible configuraion is given by 

X = {x u ...,x N }, x 1 <---<x N , (xi G Z). 

Thus the x« are the occupied sites. We assume an initial configuration Y = {yi, . . . , Un}, 
and obtain formulas for 

(1) Py(X; £), the probability that at time t the system is in configuration x|§ 

(2) Py(x m (t) = x), the probability that at time t the mth particle from the left is 
at x. 

The formula we get for the latter extends to infinite systems 

2/i < 2/2 < > +oo. 

In particular we may take Y = Z + . (This is the step initial condition.) 

For the derivation of (1) we use the Bethe Ansatz [lj to obtain a solution of 
a differential equation with boundary conditions. The derivation of (2) from (1) 
requires the proof of two combinatorial identities. The derivation we outline for this 
is from [12] and simpler than the one in [9]. 

For step initial condition we show that F(x m (t) = x) has a representation in terms 
of Fredholm determinants. This makes asymptotic analysis possible. We assume that 
q > p, so there is a drift to the left, define 7 = q — p, and obtain asymptotics as 
t 00 for 

P 0Em(£/7) > x) (fixed m and x), 
and the limits as t — > 00 of 

P (x m (th) <~t + I 1 ' 2 s t 1/2 ) (m fixed), 

P(z m (t/7) < -c^t + c^a)^ 1 / 3 ) (m= [<rt]), 

where Ci(er) and 02(0") are certain explicit constants. 

The last limit is the distribution function F 2 (s) of random matrix theory, the lim- 
iting distribution for the rescaled largest eigenvale in the Gaussian unitary ensemble. 
These asymptotics were obtained in [2] for the case of TASEP. (That F 2 should arise 
in ASEP had long been suspected. This is referred to as KPZ universality [3].) 

2 This had been known for the case N = 2 [7]. 
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II. Integral Formulas 

1. The differential equation 

The idea goes back to pQ . There is a differential equation with boundary conditions 
and an initial condition whose solution gives Py(X;t). To state it we introduce the 
new notation u(X; t) or u(X) in place of Py(X; t)Jf] 

We first consider the case N = 2, and consider du(xi,x 2 )/dt. After an exponential 
waiting time, the system could enter state x 2 } or it could leave this state. Assume 
first that X2 — Xi > 1, so that there is no interference between the two particles. 
The system could enter the state if the first particle had been at X\ — 1 (this has 
probability u{x\ — 1, x 2 )) and moved one step to the right (probability p), and three 
other analogous ways. The system could leave the state if the first particle is at Xi 
(probability u(xi, x 2 )) and moves one step to the right (probability p) or one step 
to the left (probability q), and analogously for the second particle. These give the 
equation 

^-u(xi,x 2 ) = pu(xi - l,x 2 ) + quixi + 1,22) 
at 

+pu(x 1 ,x 2 - 1) + qu(xi,x 2 + 1) - 2u(x 1 ,x 2 ). 

But if x 2 — xi = 1 then for entering the state the first particle could not have 
been one step to the right nor the second particle one step to the left, and for leaving 
the state the first particle cannot move right nor can the second particle move left. 
Therefore in this case 

^u(x 1 ,x 2 ) = pu(xi - 1,2 2 ) + qu(xi,x 2 + 1) - u(xi,x 2 ). 

We could combine these two equations into one, but then the right side would have 
nonconstant coefficients. Instead, as in pQ, we observe that if we formally subtract 
the two equations we get, when x 2 = X\ + 1, 

= pu(xi, Xi) + qu(xi + l,Xi + 1) — u(xi, X\ + 1). 

If the first equation holds for all x\ and x 2 , and this last boundary condition holds 
for all Xi, then the second equation holds when x 2 = x\ + 1. So an equation with 
nonconstant coefficients has been replaced with an equation with constant coefficients 
plus a boundary condition. 

3 The reason is that if X = {x\, . . . , xn} then Py {X; t) only makes sense when x\ < ■ ■ ■ < xn, 
but for u{X;t) there will be no such requirement. 
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This was done for N = 2, but it holds for general N. Suppose the function u(X; t), 
defined for all X = {x\, . . . , x^} G Ti N , satisfies the master equation 

±uiX;t) 

N 

= ^2\pu(...,Xi-l,...) + qu(...,Xi + l,...)- u(X)\, 

i=l 

and the boundary conditions 

u(. . .,Xi,Xi + 1,...) =pu(. . .,Xi,Xi, ...) +qu(. ..,Xi+ 1,^ + 1, ...) 
fori = 1,...,JV- 10 Suppose also that it satisfies the initial condition 

u(X; 0) = Sy(X) when x% < • • • < Xn, 
which reflects the initial configuration Y. Then 

u(X] t) = Py(X; t) when x x < • • • < xjv- 

Thus the strategy will be: (1) find a large class of solutions to the master equation; 
(2) find a subset satisfying the boundary conditions; (3) find one of these satisfying 
the initial condition. The last will be the hard (and new) part. 

2. Solutions to the master equation 

Define 

£(0=Pf 1 + ^-l- 
For any nonzero complex numbers £i, • • • , £jv, a solution of the equation is 

Since the £j are arbitrary another solution is obtained by permuting them. Thus, for 
any a in the symmetric group another solution is 



nq„ n 



e e(£i)*. 



(The second factor is symmetric in the £j, which is why we can write it as we do.) 
Since the equation is linear, any linear combination of these is a solution, as is any 



4 For N > 3 the boundary conditions arising from configurations with more than two adjacent 
particles automatically follow from the boundary conditions arising from two adjacent particles. 
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integral (over the £j) of a linear combination. Thus we arrive at the Bethe Ansatz 
solutions 

The -F CT are arbitrary functions of the and the domain of integration is arbitrary. 

3. Satisfying the boundary conditions 

We look for functions F a such that the integrand satisfies the boundary conditions 
pointwise. The ith boundary condition is satisfied pointwise when 

F ° (P + 9^(0 ^(i+i) " &(*+!)) (&(i)&(i+i)) x< II Ctf) = °" ( 2 ) 

Define Tj<r to be the permutation that differs from o by an intechange of the zth 
and (i + l)st entries. Thus, if a = (2 3 1 4) then T 2 a =(2134). Since % is 
bijective, ([2]) is unchanged if each a in the summand is replaced by TjO", and therefore 
unchanged if we add the two. Since the last two factors are unchanged upon replacing 
a by TjO", we see that a sufficient condition that 02]) satisfied is that for each a, 

Because the expression will appear so often it is convenient to define 

so the condition becomes 

F a f(€ff(i+l), &r(i)) + ^T i( T f(€a(i), £<r(i+l)) = 0- 

This is to hold for all a and all i. Since these are (n — 1) n! linear equations in the 
n\ unknowns F a , we cannot necessarily expect a solution. But there are solutions, 
and in fact it is easy to see that 

Fa(0 = SgnO U/(6r(i), &(,)) X y>(f), 

where is an arbitrary function of the £j, satisfies the equations. (In fact this is the 
general solution, since if any F a is known then all others are determined.) These F a 
in (JTJ) give a family of solutions to the master equation that satisfy the boundary 
conditions. 
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4. Satisfying the initial condition 

The initial condition is 

/ J2 F ^UChd N Z = 5y(X). 

We begin with the fact that if C is a contour enclosing zero thei]@ 
Therefore 



(3) 



Thus if id denotes the identity permutation, then the a = id summand in (J3J) will 
give the integral 5y(X) if the integration is over C N and 

i 

For this to hold we choose 

^=n/(6.&r i -n^* 1 - 

If we define 

^"w (4) 

then the solution we have chosen is 

»(*;*) = e / ^(o n%) n (^ _i ^ te)t )A. (5) 

a ^ CN i i 

It satisfies the master equation and boundary conditions, and the a = id summand 
satisfies initial condition. 

Observe that because of the poles of A a when a ^ id this will not be well-defined 
until we specify C further. 

a. TASEP 

When p = 1 we have 

^=s g u<7jj(i-^ w r w - < . 

i 

Because of this product structure the integrals of products in ((5} may be written as 
product of integrals and © becomes 



5 Unless specified otherwise, all contours are described once counterclockwise, and all contour 
integrals have a factor l/2iri. 
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u{X- 1) = J2 sgn a TT / (1 - C^' 1 e^- 1]t d£ 

i Jc 

= det Qf (i - ^-^i-^-yr 1 -!)* _ 

Schiitz [7] obtained this solution to the master equation, using Bethe Ansatz as we 
have described, and went further to show that it satisfies the boundary condition 
when the point £ = 1 is outside the contour C. So he established the formula 

P Y (X;t) = det ^ (1 -^-^-^-yr 1 -!)^ ; 
where C r denotes the circle with center zero and radius r < 1. 



b. ASEP 



In [7J Schiitz also considered ASEP and showed that when N = 2 the probabil- 
ity Py(X]t) is equal to a sum of a two-dimensional integral and a one-dimensional 
integral. In the two-dimensional integral the contours were different. It turns out 
that if one integrates over small contours only then the sum is the sum of two two- 
dimensional integrals. And this extends to general N. 

Recall that because of the poles of A a , it matters which contours C we take in (JSJ). 
When p ^ all poles of the A a will lie outside C r if r is small enough. These are the 
contours we take. 

Theorem [91 Theorem 2.1]: Suppose p ^ and assume that r is so small that all 
poles of the A a lie outside C r . Then 

p Y (x;t) = j2 f N Mt) nc (i) n (c w " i ^*)d Ar e- (6) 

a i i 

For the proof we have to show that the initial condition is satisfied. Since the 
a = id summand satisfies the initial condition, what is to be shown is that if 

•* C ? i i 

then 

E = °^ ( ? ) 

6 Thcre was an error in the proof of this in [S]. A corrected version of [H] is posted as 
larXiv:0704.2633l We outline the proof here. 
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We show that the permutations in §Ar\{i<i} can be grouped in such a way that 
the sum of the 1(a) fron each group is equal to zero. For l<ii<Afiixn-l distinct 
numbers ii, . . . , % n -\ G [1, iV — 1], define A = {ii, . . . , z n _i}, and then 

S N (A) = {a G 8> N : cr(l) = i u . . . , a(n - 1) = i n _ x , a(n) = N}. 

When n = 1 these are all permutations with er(l) = N. When n = N — 1 each Sjy (A) 
consists of a single permutation. Let B be the complement of A U {N} in [1, N]. 

We first we make the substitution 



0i<7V 

in all the integrals. The product of the powers of the in §§§ becomes 

-l(n~ x n+VN— Hi— 1 

We use the alternative representation 



Kk 

-He)>cr-i(k) 



a.- n (-^fff) m 



to see what happens when we shrink some of the ^-contours. The only poles we 
might cross come from the denominatiors in ([9]) after the substitution, and these are 
in the ^-variables when i e B. 

Lemma 1. When n = N — 1 we have 1(a) = for a G §jv(^4). 

There is a single % G B and in this case there is no pole in the ^-variable coming 
from (jHJ). Using xn > £jv-i and > in (jSJ), we see that the integrand is analytic 
at £j = 0. Therefore the integral with respect to £j is zero. 

Lemma 2. When n < N — 1 all J(cr) with a G §at(^4) are sums of lower-order 
integrals in each of which a partial product in (Q independent of a G A is replaced 
by another factor. In each integral some with « G -B is equal to another £j- with 

ie5. 

If j = max 5, we shrink the ^-contour and obtain (iV — l)-dimensional integrals 
coming from poles associated with the variables with k G B\{j}. For each such k 
we integrate with respect to £ fc the residue at this pole by shrinking the contour, and 
obtain (N — 2)-dimensional integrals having the property described in the lemma. 

Lemma 3. For each integral of Lemma 2 there is a partition of E>n (A) into pairs 
a, a' such that 1(a) + I(a') = for each pair. 
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Consider an integral in which & = £j. We pair a and a' if cr _1 (z) = o'~ (J) and 
= and <7 -1 (A;) = a'^ 1 (k) when fc 7^ j. The factor (jSJ) is clearly the 

same for both when = and A CT and y4 CT / are negatives of each other then. 

Here is why. Assume for definiteness that % < j and u~ l {i) < a~ x {j). Then the 
factor corresponding to I = i and k = j does not appear for o in ([9]) but it does 
appear for a'. This factor equals —1 when = And it is straightforward to check 
that for any k ^ i,j the product of factors involving k and either i or j is the same 
for a and o' when £j = £j. 

Now ([7]) can be shown by induction on N. When iV = 2 it folows from Lemma 1. 
Assume N > 2 and that the result holds for N — 1. For those permutations for which 
cx(iV) = iV we integrate with respect to £1, ... , (,n-i and use the induction hypothesis. 
The set of permutations with a(N) < N is the disjoint union of the various §/v(^4), 
and for these we apply Lemmas 1 and 3. 



5. The left-most particle 



The probability Py (xi(t) = x) is the sum of P(X; t) over all X for which x± = x, 
thus over all x 2 , • • • , xn statisfying x < x% < ■ ■ ■ < xn < 00. When r < 1 we may 
sum under the integral sign in dSJ), and the integrand becomes 



C<t(2)Cct(3) ' ' ' £a(N) 



(1 - £<j(2)£<t(3) • • •£ CT (AT))(1 - £cr(3) • • ■Ccr(Ar)) " ' ' (1 - £ct(AT)) 

Fortunately we have our first combinatorial identity^ 



^(2)^ct(3) ' ' ' C(7V) 



= „AT(7V-l)/2 II<j(£j ~&) 



Therefore we obtain 



7 Doron Zeilberger saw the identity when it was still a conjecture and suggested to the authors 
that an identity of I. Schur [5j Problem VII. 47] had a similar look about it and might be proved in 
a similar way. This led to the proof we present. 
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Theorem [HI Theorem 3.1]: If p ^ and r is so small that all poles of the A a lie 
outside C r , then 

= *) = ^ (JV - 1)/2 / II 4^ nn"n II(^-V^) 



Identity f[TU|) is proved by induction on iV. Call the left side of the identity 
<Pn{£i, ■ ■ ■ , £at) and the right side 4>n(£,i, ■ ■ ■ , £jv)> an d assume the identity holds for 
N — 1. We first sum over all permutations such that a(l) = k, and then sum over k. 
If we observe that the inequality i < j becomes j ^ i when i = 1, we see that what 
we get for the left side of (fT0|) . using the induction hypothesis, is 

1 - 

— — - — - V(-i) fe+i n TT & ^at-xCCi, • • • , e fc -i, • • • , 

If we substitute for ^n-i(£i, . . . , £fc-i,£jM-i> • • • , Cat) what it is and do some algebra, 
we find this would equal the right side of (ITU]) if a simpler identity held: 

N N -t N-1 

This one is proved by considering the integral 



TT(p + ^ - 2;) • -j; r • N 

ii z(p-qz) niife- 



over a large circle. The integral, and so the sum of the residues at 0, the and p/q, 
equals zero. This sum is equal to the difference of the two sides of (TTTj) . 



6. The general particle 



The probability Py(x m (t) = x) is the sum of P(Jf;i) over all X for which x m = 
x, thus over all x±, . . . , x m -\ statisfying —00 < x\ < ■ ■ ■ < x m -i < x, and all 
x m+ i, . . . ,xn satisfying x < x m+ i < ■ ■ ■ < xjy < 00. The latter we can do, as in 
the last section, since r < 1. Eventually we shall expand the ^(^-contours when 
i < m to Cr with R > 1 so that we can sum over these 

First take a partition (S-, S + ) of [1, N] with \S-\ = m — 1 and sum over all those 
a for which cx([l, m — 1]) = S_ and er([m, N]) = S+. (At the end we will sum over 
these partitions.) Set ct_ = cr|[i, TO -i]j a + = cr| [ m , jv]- Then cr_ may be associated in 
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an obvious way with a permutation in S m _i and a + with a permutation in S^-m- In 
particular, sgn a± make sense, and counting inversions shows that 

sgn a = (— l} K ( s -> s +> sgn <7_ sgna + , 

where we define in general 

K{U,V)=#{(i,j):ieU, jeV, i>j}. 

When we write 

UCh = UCh ne (i) , (12) 

the a in the first product on the right may be replaced by a_ and the a in the second 
product on the right may be replaced by a + . 

Similarly, may rewrite fll]) as 

n /(£*•)» £*■)) n n /(^e^) 

i<j<m i£S-,jeS+ m<i<j , , 

and the cr in the first product may be replaced by cr_ and the a in the last product 
may be replaced by a + 

When we sum rii>mOi) over ^ ne Xi w ^h i > m (which we may do since r < 1) 
we obtain 



?<r(m+l)S .( m+ 2) ' ' ' ^cr(iV) 
(1 - Co-(m+l)Co-(m+2) ' ' "^(AoX 1 ~ £<r(m+l) * * ^(iV)) •••(!- £o-(AT)) ^1 



We then multiply by the last factor in the numerator in (1131) times sgn cr+ and sum 
over <t + . The result is, by (ITU1) . 

( x - n 6) n (6- - w 

p (N-m)(N-m+l)/2 ^ eS + TT gc 

What remains from (|12p and (fT3"j) is 

n /(^(^u-)) n 

i<3 



3 In the cited papers we used the notation a(U, V). 
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The next step is to expand all contours to Cr with R > 1, so that we can sum 
over the Xi with i < m. When we expand the contours we encounter poles from 
the denominators in ( JTBT) and f |T4l) . and integrating the residues would give lower- 
dimensional integrals. In fact the lower- dimensional integrals from ( IT5l) are zero 
and the lower-dimensional integrals from ( |T4l) are of the same type but with fewer 
variables. Let us see why this is so. We assume p, q ^ 0. 

We first consider the poles coming from the denominator in ffT5l) and begin by 
expanding the £tv contour to a circle Cr with R very large. From the denominator in 
(01 we encounter poles at 

S-/V — — z — 

with i < N. As in the proof of Lemma 1 we find that the residue at this pole is 
analytic for inside C r , so the integral with respect to & equals zero. (An important 
point is that the term £(£i) + £ (£7v)> which appears in the exponential in the integrand, 
becomes analytic at £j = after the substitution.) After expanding the £/v contour 
to Cr we expand the ^_i-contour. Now from the denominator in f|T5|) we have poles 
at 

QN-l — Z — 

with i < N — 1. As before, the integral with respect to & of the residue is equal to 
zero. There is also the pole at 

P 



N-l 



i - q^N 

But £ N G Cr, and if R is chosen large enough this pole is inside C r and so is not 
crossed in the expansion. 

Continuing, we find that when we expand all the contours the poles of A a do 
not contribute. But the poles of (THj) do contribute, and we get a sum of lower- 
dimensional integrals, one for each subset S' + G S + . These are are minus the integrals 
of the residues at the = 1 with k G S + \S' + . If we use 1) — p(l — and 

/(l, = q — 1) we find that this residue is a constant involving powers of p and 
g times the same integrand we had before except that there are no terms involving 
the £fc with k G S + \S' + . 

Once all contours are Cr (here R > 1 should be so large that all poles of the A a lie 
inside Cr) we may sum over all Xx, . . . , x m _i statisfying — oo < xi < ■ • • < x m _i < x. 
The result of this sum is that the product Y\ i<m in (TT5]) is replaced by 

- rr 

(6r(l) - 1) (^(1)^(2) - 1) • • • (C<7(1) • • • ^a(m-l) ~ 1) 

Now we are to multiply by sgn <j_ and sum over all a_ . An identitj{§ analogous to 







Proved by interchanging p and q in (TTUj) and letting & — > 1 /£„ 
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ffTUj) tells us that the sum equals 

n & - & 

n (m-l)(m-2)/2 '■ jeS - £2 



n (i - «. 



We now put everything together. We use the notations 

r = p/q, n{U) = k(U, Z + ) = sum of the elements of U. 
The result (a special case of [121 Theorem 3]) is that when q ^ 0, 

f> Y (x m (t) = X ) = (-l) m_1 r K (^-US;)-m(m-l)/2- m fc ? fc(fc-l)/2+(m-l)(m-2)/2 

n ( 1_ II &) 11(6- -eo 



X 



n /(& & 



n fe - 6) n( i - & 

i&SL, jes; t 

Here = and the sum runs over all disjoint sets 5_ and with \S-\ = m — 1; 
indices not otherwise specified run over £_ U £+. 

There is one more step. We take a fixed set S and first sum over all partitions 
(S-, S' + ) of S with \S-\ = m — 1. (At the end we will sum over all these S). The 
only terms that involve S- and S' + individually combine as 



ii If(>-n6): 

ieS-,jeS<L Sj Sl ie5i 



where denotes the complement of S- in S. Identity (1.9) of [9] (with slightly 
different notation) is 



E II f^-rHH™" 1 

|s_|=m-i ig5_.ie5= ^ Si jeSi 
s_cs 

where the r-binomial coefficient [^] is defined by 



|5-|-1 
m — 1 



i-n& 



(16) 



"iV" 

n 



1 — T N ) (1 — T 



;i - r) (1 - r 2 ) • • • (1 - f- 
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Hence, after some algebra, the result becomes 



Theorem (9j Theorem 5.2]: We have when q ^ 0, 

F Y (x m (t) = x m ) = (-l) m_1 r" 1 ^ 1 )/ 2 J2 r R{s) - mk q k{k - l)/2 



\S\>m 



k-1 
m — 1 



n 



i lie 



(17) 



where now k = \S\ and all indices in the integrand run over S. 
Identity (IT5|) depends on a simpler identity, 



e n 

|S_|=m-l ieS-,j€St 
S_CS 



6 - 6 



m — 1 



This is proved by induction |5*|. We first observe that the left side is a polynomial 
in the £j. The reason^) is that it is symmetric in the £j, and if we multiply it by 
the Vandermonde n^te — £j) we obtain an antisymmetric polynomial which is, 
therefore, a polynomial times the Vandermonde. Since the left side it is bounded as 
each £j — > oo it is a constant. Using the induction hypothesis and a recursion formula 
for the r-binomial coefficients we see by setting some £j = 1 that the two sides of the 
identity agree. 

To prove ffl6|) by induction, we see now that the left side is a polynomial of degree 
at most one in each and, using the induction hypothesis and the recursion formula 
for the r-binomial coefficients, that the two sides agree when any £j = 1. Therefore 
the difference is of the form c Yl^i — 1). To show that c = we use (JTHjl to see that 
after dividing by some & the two sides of (1IB]) have the same limit as & — » oo. 

III. Predholm Determinant Representation for Step Initial Condition 



1. Series representation 

Until now we assumed a system of finitely many particles. Because we can take 
arbitrarily large R in (1171) the result extends to initial configurations 

2/i < 3/2 < — — > +oo, 
where the sum is taken over all finite sets S C Z + . 
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We learned this argument from Anne Schilling. 
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For step initial configuration, where Y = Z + and yi = i, we may sum over all 
sets S with \S\ = k and so obtain instead a sum over k > m. Before that, instead 
of indexing the variables in the integrand by S we index them by {1, . . . , k}, so that 
we can sum under the integral signs for all S with |,S| = k. If S = {si, . . . , s&} with 
si < ■ ■ ■ < Sk then after renaming the variables the factor HiC^ m the integrand 
becomes Yli 4T Si aric ^ becomes Yli rSl ■ These are the only terms in (fT7|) that 
involve the individual Sj. Summing the product of these two over all satisfying 
1 < Si < ■ ■ ■ < s k < oo gives, when R> r, 

' (19) 



Wt)(6/t) ■ ■ ■ (&/r) - l) ((6/r) ■ • ■ (&/r) - l) • • ■ ((&/ 
The factor Yli<j f(£u O) -1 m the integrand may be written 



i>j 



n 

*¥=5 

If we multiply (1191) by the numerator here the rest of the integrand is antisymmetric 
in the £j. Thus the integral is unchanged if this product is antisymmetrized. We make 
the substitution & — > r/£j, use identity ( TTUj) with p and g interchanged, and find that 
the antisymmetrization is 

rife - 

i 

Thus we obtain, 

Theorem (9j Corollary to Th. 5.2]: For step initial condition we have when q ^ 0, 



w*) < *) = (-!) m E * 



fc! 

k>m 



k — 1 

k — m 



T m(m-l)/2-mk+k/2 f p q\k 2 /2 



X 



Notice that on the left side we have F(x m (t) < x) rather than F(x m (t) = x) and 
on the right side the sign is different and the factor 1 — Yli 6i is gone. This is the 
result of summing the formula for F(x m (t) = x) from — oo to x. 
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For TASEP with p = only the term k = m is nonzero, the multiple integral is 
an m x m Toeplitz determinant, and we get 



F(x m (t) < x) = det (j ^-i+v- 1 (f - \)-' m e { t- l)t d£J . 

This was obtained by Rakos and Schiitz [6] who showed it was equivalent to Johans- 
son's result mentioned in the introduction. 

2. Fredholm determinant representation 

If we make the change of variables 



then 



n 



- 6 



1 - r?7j 

1 - Vi 



'1 + r) fc ^ -1 ^ TT — — . 



The right side can be represented in terms of the Cauchy determinant 



det 

Going back to the & gives the identity 

6 - 6 



r 



fc(fc-i)/ 



1)/2 n^-fa-^) 



n 



p + - & 



1)* ( pg )-^-D/2 JI(l-6)(g&-p)-det 



1 



p + g&£» - 6 



J S*/ l<i,j<fc 



The theorem becomes 

P(x m (t) < x) 



k>m 



k-l 

k — in 



T 



(l-m)fc 



X ■ 



det(K(e 4 , &)) l<i j'<fc 



where 



i p e(£)t 



p + ?ec - £ 

Denote by K the operator acting on functions on Cr by 

Jc R 
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The Fredholm expansion is 



det(I -\K) = J2 
which gives 



-A)' 



k=0 



k\ 



det K(£i, ij)i<ij<k d£i-~ d£ k 



k\ 



det K(£i, £j)i<i,j<k d£i--- d£ k 



where we take any contour enclosing A = 0. Thus, 



F{x m {t) <x) = (-l) m T m(m - 1)/2 

fc>m 



det(J - XK) 
W 1 



dX, 



' k 


- 1 ' 


T (l-m)k 1 


_k 


— m 


T J 



det (J - XK) 
W 1 



dX. 



If the contour is C p with p > r 1 m then we can interchange the sum and integral and 
use the r-binomial theorem 



k>m L 



k ~ 1 

k — m 



n r 



with 2 = t m A . We obtain, 



Theorem (TOj Formula (1)]: We have when q ^ 0, 

m . , /• det(I-XK) dX 
F(x m (t) < x = / ^- — , 

where the contour of integration encloses all the singularities of the integrand. 



(20) 



We can evaluate the integral by residues, getting a finite sum of determinants. 
When m = 1 we obtain 

P(zi(£) > x) = deb(I-K), 

IV. Asymptotics 

1. Statements of the rsults 

If q > p and t — > oo, we expect x m (t) to be large and negative. We obtain three 
asymptotic results for P (x m (t) < x) as t — > oo. Recall the definition 7 = q — p. 

Theorem 1 [TU Theorem 1]: Let m and x be fixed with x < m. Then as t — > 00 

00 

P(x m (t/ 7 )>a;)~ J](l-r A 



-j-2m—x—2 g— 4 



fe=l 



(m — 1)! (m — x — 1)! 
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Theorem 2 [TTJ Theorem 2]: For fixed m we have 
where i^, is the operator on L 2 (s, oo) with kernel 

z') = q e -(P 2 +q 2 ) (z 2 +z' 2 )/4+pqzz' 

For the third result we recall that 

F 2 (s) = det (7 - K Aily X( s ,oo)), 

where 

roc 

K A iry{x,y)= Ai(x + z) Ai(y + z) dz. 
Jo 

For a G (0, 1) we set 

ci(a) = -1 + 2v^, c 2 (a) = a" 1 / 6 (1 - v^) 2 / 3 . (21) 

Theorem 3 (TTJ Theorem 3]: We have 

lim P (x [at] (t/j) < ci(o-) t + c 2 (a) st 1/3 ) = F 2 (s) 

uniformly for o in compact subsets of (0, 1). 

For TASEP (p = 0) the probabilities are m x m determinants. For m and x fixed 
the asymptotics of the determinant are easily found and agree with Theorem 1. 

A special case of Theorem 2 is 

limP(xi(t/ 7 ) > -t-7 1/2 st 1/2 ) =det(I-K B ). 

t—too 

This is a family of distribition functions parameterized by p G [0, 1). When p = it 
is a normal distribution and the probability on the left is the probability for a free 
particle. 

Theorem 3 when p = gives the asymptotics for TASEP obtained by Johans- 
son [2]. A consequence of Theorem 3 is that for fixed s G (0, 1) 

r x [<?t]{t/l) , , 
hm — — = Ci(er). 

t— >oo t 

in probability. In fact Liggett [I] showed that this holds almost surely. 
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2. Preliminaries 



A natural approach to the asymptotics is to look for a limiting operator such 
that det(J — XK) — > det(J — XK^) as rj — >■ oo. Once one has guessed i^oo there are 
two possible approaches: 

(i) Show that K — > in trace norm. 

(ii) Show that tr/T 1 — > tr for each n G Z + and that K is bounded in trace 
norm (or even Hilbert-Schmidt norm) as t — > oo. This suffices because of the general 
formula 

logdet(J — XL) = — ^ — trL™, (22) 

n=l n 

which holds for sufficiently small A. 

Both approaches will be used eventually. The operators K on Cr have exponen- 
tially large norms as t — > 00, and we will replace them by operators with the same 
Fredholm determinants that are better-behaved. This will be possible because of 
lemmas on stability of Fredholm determinants. 

Lemma 1. If s — > T s is a deformation of closed curves and L(rj, rj') is analytic for 
77, 7/ G r s for all s, then the Fredholm determinant of L acting on T s is independent 
of s. 

Lemma 2. If Li(rj, rj') and L 2 (rj, rj') are two kernels acting on a simple closed 
curve T, such that Li(r], ?/) extends analytically to rj inside V or to rj' inside T, and 
-^2(7/) v') extends analytically to r] inside V and to v( inside T, then the Fredholm 
determinants of Li(r), rj') + L 2 {ri, rj') and la (77, rj') are equal. 

The proofs use the fact that det(J — XL) is determined by the traces trZA For 
Lemma 1 we use 

trL n = ■■ L(rji, rj 2 ) ■ ■ ■ L(r] n -i, i] n ) L(r] n , 771) drji ■ --dr/n. 

By analyticity the integral is invariant under the deformation. For Lemma 2, we have 
to show 

tr {L l + L 2 ) n = ti VI. 
If, say, Li (rj, rj') extends analytically to 77' inside T, then 

L 1 L 2 (rj, rj") = J L,(rj, rj') L 2 (r,' , rj") dr,' = 0. 

Similarly L\ = 0. Also, trL 2 = 0, so tr (L x + L 2 ) = trl^. Since L X L 2 = L\ = 0, we 
have for n > 1 

(Li + L 2 ) n = LI + L 2 LT 1 . 
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Since 

trLsL™- 1 = trL™" 1 ^ = 0, 
we have tr (Li + L 2 ) n = tr L™. 

3. Another operator 

We introduce the notation 

v>{v) = 



1-tt]\ r_t i_] 

1 — 7] 



In K(£, £') we make the substitutions 

1 — rri 1 — rrf . 

£ = -; £' = -; r> *^*/7 

1 — 77 1 — 7]' 

and obtain the kernel"] 

77' — Try 

acting on c, a little circle about 77 = 1 described clockwise, which has the same 
Fredholm determinant. We denote this by K 2 because there is an equally important 
kernel 

7] — T7] 

The kernel ifi(ry, ry') extends analytically to 7] and ry' inside c while K 2 (r], rf) 
extends analytically to 77 inside c. Hence by Lemma 2 the determinant of K 2 equals 
the determinant of K 2 — K\ . 

Next we apply Lemma 1 to the kernel 



7/ — TT) 

with r = — c and T\ = C p with 1 < p v < r _1 . (Recall that c was described clock- 
wise.) Since the numerator vanishes when the denominator does, the only singular- 
ities of the kernel are at 7/, 7/ = 1, r _1 , neither of which is passed in a deformation 
r s , s G [0, 1]. Therefore the operator K acting on Cr may be replaced by K\ — K 2 
acting on C Pv . 



^his is the kernel {d^/drf) 1 / 2 ^' /dr]') 1 ' 2 K{^{rf), f (r/)) 
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4. Another Fredholm determinant representation 



The function (p{rrj) is analytic on sC p when < s < 1. The denominator rj' — ttj 
in Ki is nonzero for 77, 7/ e sC p for all such s. Therefore by Lemma 1 the Fredholm 
determinant of K\ on C p is the same as on sC p . This in turn is the same as the 
Fredholm determinant of 

sK^sr], ST] ) = (23) 

7]' — T1] 

on C p . As s — > this converges in trace norm to the kernel 

K (V, rf) 1 



Tj 1 — TT] 



on C p . Therefore the Fredholm determinant of K\ equals the Fredholm determinant 
of A Q . 

The kernel of Aq equals 



K 2 (v, rf) 



(C - rr]) (rf - t() r]'-T 2 T] 



because ttj is inside C p and r ^ outside when 77, 77' e C p . In particular tr Aq = 
(1 - r 2 )- 1 . Generally, we find that tr A™ = (1 - r 2 )~ n . Thus by Q22) we have for 
small A 

iogdet ( / - AA' o) = - y. -rh = - E E ^ = E - *A 

n=l k=0 n=l k=0 

and so 

00 

det(I - AAi) = det(/ - AA ) = - Ar fc ). 

k=0 

We factor out / — XKi in 

det(J-AA) dX f det(J - AA X + AA 2 ) dA 



(recall the substitution t — )■ t/7) and obtain 

/ OO 7 -v 

J] (1 - A r fc ) det(J + AA 2 (J + A)) — , (24) 

fc=m 

where A is the resolvent operator AAi (I — AAi) -1 . 
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5. Theorems 1 and 2 



Consecutive integration shows that for small A the resolvent kernel has the nice 
representation 

n=l ' ' 

where 

<fn(v) = V>{v) <f( T v) ■ ■ ■ ( f(T n - 1 7]). 

For Theorems 1 and 2, whose derivations we shall not explain in detail, we wrote R = 
R 1 +R 2 where R\ is analytic everywhere except for the poles at A = 1, r _1 , . . . , r~ m+1 
and R 2 is analytic for |A| < r~ m . For Theorem 1 the asymptotics comes from the 
residue of Ri at A = r - m +! p or Theorem 2 we used approach (ii) described above. 
In [10] a steepest descent computation had shown that tr K n — > tr K n for all n. What 
was needed to complete the proof was to show that K 2 (I + R) had bounded Hilbert- 
Schmidt norm as t — > oo, uniformly for A in compact sets not containing any of the 
singularities r~ k . We used the representation R = Ri + R 2 to show that this was so. 



6. Theorem 3 



Here m = at is large and 1, r 1 , . . . ,r m+1 must be inside the contour. If we set 
A = \i r~ m we can take ii G C p with p > r fixed, and (I2TI) becomes 

P(x m (t/ 7 ) < x) = [ f[(l - n r k ) det(J + li r~ m K 2 (I + R))^L. (26) 
In f )25|) we use 



where 

V?oo(?7) = lim cp n (r)) = (1 - i)Y x e~ l . 

n— ¥oo 

The Cauchy integral representation of y?oo(T n rj)~ l , and some manipulation of series 
and integrals, give 

K 2 (I + R) ( V , rf) = -( 77 ^9, ^ rf C 



OO 



Op. (C - "7) (V - C) 



' 1 — r ft A /r C — TV 



j.. _. x - r ^ ic Pc C - TV J Cpu Voo(u) (rf - u/t) U k+l ' 
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where the radii of the contours in the series satisfy 

The first operator on the right side is analytic for |^|, |r/'| < p v and the infinite 
sum is analytic for \rj\ < p v . It follows by Lemma 2 that the Fredholm determinant 
of the sum of the two, i.e., of K%{I + R), equals the Fredholm determinant of the 
infinite sum. 

If we set 



T 



k 



f(ji,z)= 7 



..k 



k=—oo 



then since A = pr m , 

-k / /-\ k 



k=—oo 

and so the infinite sum may be written 



^(;) =^g)"/o,c/«). 



c Pu Je P( Vqo(m) W (( ~ ti]) (v' ~ u/t) u 
The substitutions 77 — )■ 77/1", 77' — > 7//t replace this by the kernel 

^oo(C) fCV f(j*,(M du 



c Pu Jc n <Poo(u) \uj (C-v)(v'-u) u 

where now the operator acts on with 6 (r, 1) and in the integral 

p c e (1, r -1 ^!;), P« e (rp c , p,). 

If we expand the u-integral so that p v < \u\ < 1 on the new contour we pass the 
pole at u — which gives the contribution 

vUC) C m f(^C/v') dc (27) 



PC 



<pM) (v') m+1 C-v 



The new double integral is analytic for \r)\, \r)'\ < p v and (l27|) is analytic for |?7| < p v . 
Therefore by Lemma 2 the Fredholm determinant is the same as that of (I27I) . 

We have shown that if we define 

•/c„ c <Poo(r/') (77') m+1 C - 7? 
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where p( 6 (1, r 1 Pv), then fl26|) becomes 

P(x m (t/ 7 ) <x)= I J~J(1 - /J r k ) det(I + /i J) — . (29) 

fc=0 ^ 

This representation, in which the parameter m is in the operator, makes an asymptotic 
analysis possible. 

By Lemma 1 the contours C Pi] (the home of the functions on which J acts) and C Pc 
(in the integral defining J) may be simultaneously deformed if during the deformation 
we do not pass a singularity of the integrand. 

We apply steepest descent, and so look for the saddle points for </?(£) ( m when 
m ~ at and x ~ ct. In general there are two saddle points. When c equals c\(a), 
given in (l2~Tj) . they coincide at 

£ = -V^/(l-V^). 

Both contours may be deformed to pass through the saddle point, the neighborhood 
of which gives the main contributions. If x = C\(a) t + 02(a) st 1 ^ 3 precisely, and we 
make the variable changes 

for a certain constant C3, then the rescaled kernel /i J(p, //) has limit 

e -C 3 /3+<+(V) 3 /3-sr,' 

e (C - 17) (1/ - dC ' 

(The constants 02(a) and C3 come from a third derivative at the saddle point.) Here 
I\ consists of the the rays from to 00 e ±27rl//3 . The limiting operator acts on functions 
on the contour Y v consisting of the the rays from to 00 e ±7T ^ 3 . 

Using the general identity det(J — AB) = det(I — BA) we replace this by the 
kernel 



e -C 3 /3+»? 3 /3+yC-^ 



C-77 



dr/d( = -K AiTy (x, y), 



acting on L (s, 00), where 



/"OO 

^Air y (a;, y)= Ai(z + x) Ai(z + y) dz^ 
Jo 



Hence 

det(J + n J) det (J - K Airy X( S)0 o)) = i^aOO 
for all /x, and it follows that the integral in (1291) has the limit ^(s) 



12 The reason the double integral equals —if Airy (2;, 2/) is that applying the operator 9/9cc + d/dy 
to the two kernels gives the same result, Ai(x) Ai(y), so they differ by a function of x — y. Since 
both kernels go to zero as x and y go to +00 independently this function must be zero. 
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